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ABSTRACT 

The mass of massive black holes in quasar cores can be deduced using the typical 
velocities of H/3-emitting clouds in the Broad Line Region (BLR) and the size of this 
region. However, this estimate depends on various assumptions and is susceptible to large 
systematic errors. The H/3-deduced black hole mass in a sample of 14 bright quasars 
is found here to correlate with the quasar host galaxy luminosity, as determined with 
the Hubble Space Telescope (HST). This correlation is similar to the black hole mass vs. 
bulge luminosity correlation found by Magorrian et al. in a sample of 32 nearby normal 
galaxies. The similarity of the two correlations is remarkable since the two samples 
involve apparently different types of objects and since the black hole mass estimates in 
quasars and in nearby galaxies are based on very different methods. 

This similarity provides a "calibration" of the H/3-deduced black hole mass estimate, 
suggesting it is accurate to ±0.5 on log scale. The similarity of the two correlations also 
suggests that quasars reside in otherwise normal galaxies, and that the luminosity of 
quasar hosts can be estimated to ±0.5 mag based on the quasar continuum luminosity 
and the H/3 line width. Future imaging observations of additional broad-line active 
galaxies with the HST are required in order to explore the extent, slope, and scatter of 
the black hole mass vs. host bulge luminosity correlation in active galaxies. 

Subject headings: galaxies: nuclei-quasars: general 
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1. INTRODUCTION 

Indirect evidence for the existence of massive black 
holes (MBHs) in Active Galactic Nuclei (AGNs) has 
been growing over the years (e.g. Rees 1984). How- 
ever, the most conclusive evidence for the existence 
of massive black holes has been recently obtained in 
the Milky Way (e.g. Gcnzcl ct al. 1997), and in 
NGC 4258 (Miyoshi et al. 1995), a weakly-active 
galaxy. This new evidence is based on high spatial 
resolution observations of stellar and gas kinematics. 
Similar estimates could not be employed in quasars 
and bright Seyfert galaxies as the stellar kinematics 
close to the black hole is hopelessly lost behind the 
glare of the active nucleus. A rough estimate of the 
black hole mass in AGNs can be obtained based on 
the size and the typical velocities in the Broad Line 
Region (BLR, e.g. Dibai 1981; Wandel & Yahil 1985; 
Joly et al. 1985; Padovani & Rafanelli 1988; Koratkar 
& Gaskell 1991). However, this method is susceptible 
to various systematic errors, and there is currently no 
independent way to estimate its accuracy. 

Compact massive dark objects, most likely MBHs, 
were inferred in the cores of many nearby normal 
galaxies based on stellar kinematics and the observed 
light distribution (see review by Kormendy & Rich- 
stone 1995). In a recent comprehensive study of the 
stellar dynamics of a large sample of nearby galax- 
ies Magorrian et al. (1998) found that a MBH may 
exist in the cores of nearly all bulges. They also con- 
firmed the strong correlation between the black hole 
mass and the bulge mass, consistent with Mbh ~ 
0.006Mbuigc- If quasars reside in normal galaxies, 
then their black hole mass and bulge mass should 
also follow this correlation. It is not possible to ex- 
plore this correlation directly in quasars since stellar 
velocity distributions in the host bulges have not been 
measured yet. However, since there is a strong corre- 
lation between Mb u i ge and Lbuigc in galaxies (Faber 
et al. 1997), one can instead test if quasars follow 
the Mbh versus Lbuige correlation found for normal 
galaxies. 

Rather accurate determinations of quasar host galaxy 
luminosities were recently obtained by Bahcall et 
al. (1997) for a representative sample of 20 bright 
low rcdshift quasars using the Hubble Space Tele- 
scope (HST). In this Letter I show that the quasar 
host galaxy luminosity appears to be significantly 
correlated with the H/3-deduced black hole mass, 
AfBH~(H/3), and that this correlation is very similar to 



the ibuige versus Mbh relation determined by Magor- 
rian et al. for nearby normal galaxies. This similar- 
ity provides a "calibration" for the M B h estimates in 
AGNs. The Mbh(H/3) estimation method, its appli- 
cation to the Bahcall et al. sample, and the correla- 
tion analysis results are given in §2, and the impli- 
cations are discussed in §3, together with some pre- 
dictions of host luminosities which can be tested with 
HST in the near future. 

2. The Mbh(H/3) versus Lh os t Correlation 

2.1. Af BH (H/3) 

The black hole mass can be estimated using the 
the velocity dispersion of the H/3 emitting clouds in 
the BLR and the size of this region, together with the 
assumption that the clouds' motion are virialized, i.e., 

M BH (H/3) = i?BL R (H/3)4 LR /G (1) 

where -Rblr(H/3) is the size of the H/3-emitting region 
in the BLR and wblr is the observed H/3 velocity 
dispersion. 

Kaspi et al. (1996) find J? B lr(H/3) = O.OWi^ 2 pc, 
where L44 is the 0.1 — \[ira luminosity in units of 
10 44 erg s _1 , assuming Q = 1.0, H = 75 km s _1 Mpc -1 . 
This relation is equivalent to 

i?BL R (H/3) = 0.086L^ 2 pc, (2) 

where L4Q is the Bolometric luminosity in units of 
10 46 erg s" 1 (using L Bo i = 3L .i-i M m, e.g. Fig.7 in 
Laor & Draine 1993), and Q = 1-0, H = 80 km s" 1 Mpc" 1 
which is used throughout the paper. The -Rblr oc 
L 1 / 2 scaling is also indicated by the weak luminosity 
dependence of AGN emission line spectra, and is also 
expected based on dust sublimation which occurs at 
i? dust ~ 0.2L^ 2 pc (Laor & Draine 1993; Netzer & 
Laor 1993). 

Thus, just based on vblr, taken here as the ob- 
served H/3 FWHM, and L46, one obtains the following 
estimate for the black hole mass from Eqs.l& 2, 

m 9 = 0.18A4, 00 i^ 2 , (3) 

where m 9 = M BH (H/3)/lO 9 M , and 
AW3000 = H/3 FWHM/3000 km - 1 . 

2.2. Lhost 

Quasar hosts have been studied extensively from 
the ground (e.g. McLcod & Ricke 1994a, 1994b; Dun- 
lop et al. 1993). However, separating out the quasar 
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host galaxy clearly requires a high angular resolution, 
and thus measurements with the iJSTcan provide the 
most accurate determination of quasar host proper- 
ties. 

I use the sample of 20 luminous low redshift (z < 
0.3) quasars studied by Bahcall ct al. with the HST 
Wide Field/Planetary Camera-2 (WFPC2). This 
sample is likely to represent the properties of nearby 
bright quasars. In addition, the large sample size, the 
uniform and detailed reduction, and the detection of 
all quasar hosts, make this sample the best one avail- 
able for exploring the Mbh versus Lhost relation in 
quasars. 

Some of the host galaxies morphologies were iden- 
tified by Bahcall et al. as elliptical, and for these 
galaxies Lbuigo = -^host (taken from the best fit 2- 
D model in their Table 5). Other hosts were iden- 
tified as spirals, or interacting, and the value of 
-^buigo for these objects, required for a direct com- 
parison with the Magorrian et al. results, is not 
available. An estimate of Lbuigo for the objects 
which arc best fit by an exponential disk is ob- 
tained by subtracting the 7.5 < r < 15 kpc annular 
magnitude (their Table 8) from the total magnitude 
(their Table 5), yielding My (inner host). The mean 
AMy(host— inner host) is 0.5 mag, which is smaller 
than (AMb (total— bulge)) ^1 — 2 mag for early type 
spiral galaxies (Simien & de Vaucouleurs 1986), sug- 
gesting that My (inner host) overestimates Lbuige- Al- 
though, My (host) may be underestimated, since the 
fit does not included a bulge component. 

2.3. The correlations 

Table 1 lists the Bahcall et al. quasars used for the 
correlation analysis together with their z, My (bulge) 
for objects with a de Vaucouleurs fit, or My(inner 
host) for objects with an exponential disk fit, H/3 
FWHM, bolometric luminosity, Mbh (H/3) as deduced 
from Equation 3, and the host morphology from Bah- 
call et al. The H/3 FWHM is obtained from Boroson 
& Green (1992) which provide high quality and uni- 
formly reduced spectra of all 87 z < 0.5 PG quasars 
(Schmidt & Green 1983), of which 14 overlap with the 
Bahcall et al. sample. Continuum fluxes are avail- 
able for all of these 14 PG quasars in Neugebauer 
et al. (1987), which provides accurate and uniformly 
reduced continuum spectrophotometry for most PG 
quasars. The luminosity at 3000A is converted to 
LboI using L Bo i = 8.3 x iaL„(3000A) (see Fig.7 in 
Laor & Draine). 



There is no uniform data set with the continuum 
flux and H/3 FWHM for 5 additional quasars from 
the Bahcall et al. sample. Different papers quote pa- 
rameters that can differ by > 50% for a given object. 
These objects were therefore not included in the anal- 
ysis as they may be subject to significant systematic 
deviations. 

The upper panel in Figure 1 shows My (bulge/inner 
host) [hereafter My (b/ih)] versus Mbh for the 19 
Bahcall ct al. quasars. Only the 14 quasars marked 
with filled squares were used in the analysis. The 
Spearman rank order correlation coefficient is —0.70 
which has a probability of 0.005 to occur for unre- 
lated parameters. A simple least squares fit to the 
data gives 

My(b/ih) = -21.76±0.24-(1.41±0.38)logm 9 . (4) 

Three quasars which are best fit by an exponential 
disk, PKS 1302-102, PG 1307+085, and PG 1444+407 
(Bahcall et al. Table 5), may have an elliptical mor- 
phology (Table 1). A least squares fit for the 14 
quasars using the three quasars de Vaucouleurs fit 
My (bulge) yields the coefficients -21.85 ± 0.28 and 
-1.18 + 0.44. 

The middle panel in Fig.l shows the My (bulge) 
vs. Mbh relation for nearby normal galaxies from 
Magorrian et al. (1998). The dashed line represents 
the relation 

My(bulge) = -21.40 - (2.21 ± 0.28) logm 9 , (5) 
as deduced from the 

log(M BH /M ) = -1.79+(0.96±0.12)log(M bulgo /M Q ), 

(6) 

and 

log(M bulgc /M ) - -1. 11+(1. 18+0.03) log(L bulgc /L ) 

(7) 

relations found by Magorrian et al. , and the standard 
relation My(bulge) = 4.83 — 2.51og(Lbuige/io)- 

The quasar correlation is flatter than the Magor- 
rian et al. correlation (—1.41 vs. —2.21). This may 
be partly due to the fact that all the quasar hosts 
at logTOg < —0.7 are disk galaxies, where My (inner 
host) may overestimate My (bulge) (§2.2). 

The lower panel in Figure 1 compares directly the 
distributions of the 32 Magorrian et al. galaxies and 
of the 19 quasars in the Mbh versus My (b/ih) plane. 
The two distributions overlap surprisingly well. 
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3. DISCUSSION 

The overlap of the distributions of quasars and of 
normal galaxies in the Mbh versus My(b/ih) plane 
is the main result of this paper. This overlap is re- 
markable as bright quasars and nearby galaxies are 
apparently different types of objects, and since the 
Mbh estimates in quasars and in nearby galaxies are 
based on very different methods (the BLR versus stel- 
lar dynamics). The overlap is also surprising given 
the crudeness of the Mbh estimates for both pop- 
ulations. As stressed by Magorrian et al., their data 
are fit with a simplified, axisymmctric, stellar dynam- 
ics model, and a more general model may yield Mbh 
which could be off by an order of magnitude, or may 
even not require a "massive dark object" at all. The 
Mbh (H/3) estimate is even cruder. Large systematic 
errors could be induced if the BLR velocity field or the 
optical-UV continuum are anisotropic, if the scaling 
of i?BLR.(H/3) with L does not hold in bright quasars, 
or if the H/3 dynamics is affected by non gravitational 
forces (e.g. radiation pressure, magnetic fields). 

The overlap suggests a number of interesting im- 
plications. First, concerning the BLR: 1. The H/3 dy- 
namics are most likely dominated by gravity; 2. the 
H/3 velocity field and the observed optical-UV emis- 
sion are not likely to be strongly anisotropic, and 3. 
the -Rblr(H/3) versus L relation most likely holds in 
quasars. Second, concerning the Mbh (H/3) estimate; 
the overlap allows a "calibration" of this mass esti- 
mate and suggests it is probably accurate to within 
±0.5 on log scale. Third, concerning quasar hosts; 
1. the scatter in the Mbh (H/3) versus My(b/ih) cor- 
relation suggests that My(b/ih) can be estimated to 
within ±0.5 mag based on the quasar luminosity and 
H/3 line width. 2. The overlap of the two distribu- 
tions suggests that quasar hosts are similar to normal, 
nearby galaxies, and thus that My(b/ih) is generally 
not strongly affected by processes such as a nuclear 
star burst, or distortions due to a tidal interaction. 

This correlation may also be useful for surveys of 
the large scale structure of the universe (2dF, SDSS). 
Quasars can be used as bright markers of galaxies out 
to high z whose bulge luminosity and mass can be 
deduced from the quasar emission spectra, allowing 
studies of clustering as a function of mass. 

The correlation found here, M v (b/ih) cx M B ^ A±aA 
(or M B ^ 2±0 ' 4 ) translates using Eq.7 to M BH oc M*£f c 
(or M^ 8 * ' 6 ), which is steeper than the Magorrian et 
al. relation Mbh oc Mb u i go - The slope of the quasar 



relation has a relatively large uncertainty due to the 
small range in Mbh available (—1.16 < logm 9 < 
0.17), but it is interesting to note that at the high 
mass end (logmg > 0.2) the Magorrian et al. galax- 
ies appear to follow the quasar relation quite well (see 
Fig.l). At the low black hole mass end one has the 
two best Mbh estimates available, in the Galaxy and 
in NGC 4258, where logm 9 = -2.59; -1.44 (Miyoshi 
et al. 1995; Genzel et al. 1997) and My (bulge) = 
-18.4; -19.13 (Bahcall & Soneira 1980; RC2 cata- 
logue + Simicn & de Vaucouleurs 1986). These galax- 
ies follow the quasar relation significantly better than 
the Magorrian relation (see Fig.l). Thus, the data 
in the range —2.59 < logmg < 1.2 appears to agree 
better with the steeper quasar relation. The quasar 
relation is also interestingly close to the Hachnelt, 
Natarajan & Rees (1998) prediction of Mbh oc M^j?. 

There are some objects, such as M 32, which ap- 
pear to agree better with the nearby galaxies rela- 
tion (Fig.l). A number of Seyfert 1 galaxies with 
Mbh ~ 10 8 — 10 9 M Q , as deduced by reverberation 
mappings (Peterson et al. 1998), are 1-2 mag brighter 
than expected based on the quasar relation (Ho 1998). 
Subtraction of the AGN light from the host light 
would bring them closer to the quasar relation. 

The distribution of quasars in the absolute quasar 
B band magnitude Ms (quasar) versus the absolute 
host H band magnitude M#(host) plane appears to 
be bounded such that Ms (quasar) < Mh (host) (e.g. 
McLeod & Rieke 1995, their figure 6). McLeod (1998) 
suggested that the reason for this bound is that ob- 
jects where Mb (quasar) — M//(host) "have a maxi- 
mum allowed black hole mass for their galaxy mass 
and that the black hole is accreting at the Edding- 
ton rate." This idea is broadly consistent with the 
correlation found here. For example, a quasar with 
M^(host) = -25 mag typically has My(b/ih) ~ 
—21.3 mag (using (V — H) = 3.7 mag for our 14 PG 
quasar hosts, with M#(host) from McLeod & Rieke 
1994b). Equation 4 then gives logmg ~ —0.33 (or 
-0.47). A magnitude of M B (quasar) = My (host) = 
—25 mag translates to log i/L„(4400A)=45.55, and 
logLsoi — 46.5, which corresponds to 0.44 or 0.75 
of iEddington- The above estimates are rather rough 
since the M#(host) versus My(b/ih) correlation has 
a significant scatter. 

How can the analysis presented here be improved? 
The crude "inner host" estimate for the bulge lumi- 
nosity, used here for disk galaxies, can be improved by 
fitting a disk±bulge model to the HST images. This 
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may be feasible for early type spiral hosts where the 
typical effective radius of the bulge is r e ~ 1.4 kpc 
(Simien & de Vaucouleurs), or ~ 0".5 for the Bahcall 
et al. z ~ 0.2 quasars. However, it will not be fea- 
sible for late type spiral hosts, where r e ~ 0.3 kpc, 
and a sample of lower z AGNs will be required. One 
also needs to measure the H/3 FWHM and luminos- 
ity simultaneously to guard against variability, and to 
obtain a more accurate estimate of the ionizing lumi- 
nosity to use in the -Rblr(H/3) relation (Eq.2). Ob- 
scuration effects are well established in AGNs, and 
these may increase the scatter, if not accounted for. 
In particular, objects with a very narrow H/3 line may 
have their BLR partly obscured, or may be strongly 
dominated by emission from the narrow line region. 
Using the variable H/3 component profile can over- 
come such biases. 

Future observations with HST can address the fol- 
lowing questions: 1. Does the black hole mass vs. 
host bulge luminosity correlation extends to quasars 
with higher and lower black hole masses? 2. Can the 
scatter in the correlation be reduced with a more care- 
ful analysis, or is it intrinsically large, as suggested for 
galaxies? and 3. Does M B h oc M^^' s , as suggested 
here, or is M B h °c M bu i ge as suggested by Magorrian 
et al.? 

The PG quasars sample may be particularly use- 
ful for such future explorations with HST since a 
high S/N homogeneous spectroscopic data base is al- 
ready available from Neugebauer et al. and Boro- 
son & Green. Using this data set and Eqs.3 & 4 
one can predict that of the 87 z < 0.5 PG quasars, 
some of the lowest luminousity hosts should be found 
in PG 1244+026, PG 1404+226, and PG 1448+273 
(predicted My(b/ih)=-18.4 to -19.3), while some 
of the highest luminousity hosts should be found 
in PG 1704+608, PG 1425+267, and PG 2308+098 
(—22.0 to —22.3). One can also predict that the hosts 
of PG 2304+042 and PG 2209+184 should be 2-3 
magnitudes brighter than the hosts of PG 1244+026 
and PG 1448+273 respectively, although the former 
and later quasars have, respectively, similar luminosi- 
ties. 

This work was supported by the fund for the pro- 
motion of research at the Technion. Many thanks 
to Avi Loeb, Dani Maoz and Don Schneider for very 
helpful comments. 
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Table 1: QUASAR SAMPLE 



Object 


2 


My 




T c 






PG 0052+251 


0.155 


-20.89 


5.20 


46.02 


8.74 


Sb 


PHL 909' 


0.171 


-21.48 


11.0 


46.37 


9.57 


E4 


NAB 0205+02 7 


0.155 


-19.33 


1.05 


46.41 


7.55 


SO? 


PH 0923+201 


0.190 


-21.48 


7.61 


46.06 


9.09 


El 


PG 0953+414 


0.239 


-20.29 


3.13 


46.39 


8.49 


? 


PKS 1004+130 


0.240 


-22.48 


6.30 


46.39 


9.10 


E2 


PG 1012+008 


0.185 


-19.91 


2.64 


45.85 


8.07 


Int. 


HE 1029-140 / 


0.086 


-20.98 


7.50 


46.49 


9.30 


El 


PG 1116+215 


0.177 


-21.88 


2.92 


46.38 


8.42 


E2 


PG 1202+281 


0.165 


-20.98 


5.05 


45.51 


8.46 


El 


3C 273 


0.158 


-22.58 


3.52 


46.96 


8.87 


E4 


PKS 1302-102 


0.286 


-21.15 


3.40 


46.77 


8.75 


E4? 


PG 1307+085 


0.155 


-20.51 


2.36 


45.99 


8.04 


El? 


PG 1309+355 


0.184 


-21.23 


2.94 


45.83 


8.15 


Sab 


PG 1402+261 


0.164 


-19.95 


1.91 


45.96 


7.84 


SBb 


PG 1444+407 


0.267 


-20.49 


2.48 


46.13 


8.16 


El? 


3C 323.1 


0.266 


-21.48 


7.03 


46.34 


9.17 


E3? 


PKS 2135-147 7 


0.200 


-21.58 


5.50 


46.61 


9.09 


El 


PKS 2349-014 / 


0.173 


-22.58 


5.50 


46.57 


9.07 


Int. 



a Host absolute magnitude from Bahcall et al. (1997, Ta- 
ble 5), calculated for fi = 1.0, H = 80 km s _1 Mpc" 1 . 
b H/3 FWHM in units of 1000 km s" 1 from Boroson & 
Green (1992). 

c Log Bolometric luminosity in erg s _1 , based on /„ at 

rest frame 3000A from Neugebauer et al. (1987). 

d Log of black hole mass in units of Mq (see Eq.3). 

e Host morphology from Bahcall et al. 

f Object not included in the correlation analysis (see 

text). 
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Fig. 1. — Comparison of the correlations for quasars and for nearby galaxies. Upper panel: the My (bulge/inner 
host) versus Mbh correlation for quasars. The solid line is a least squares fit. Open squares represent objects which 
were not included in the fit (see text). Middle panel: the My (bulge) versus Mbh relation obtained by Magorrian ct 
al. for nearby normal galaxies. Lower panel: the two data sets overlaid. The two distributions overlap surprisingly 
well. The quasar relation is also consistent with the nearby galaxies distribution at log nig > 0.2, and the positions 
of the Milky Way and NGC 4258. 



7 



